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Numerical  forecasting  has  resulted  from 
accuracy  and  objectivity  in  weather 

“so'beL  farthered  by  Joint  work  ^atween  meteorology- 
especially  the  theoretical  dynamics  of  wea- 

weather  forecasting.  It  is  one  of  the  two  legs  of  wea 

ther  fcrecastiiigp 

The  demand  for  weather  forecasts 
and  the  objectives  of  forecasting  have  increased  in  number. 

Numerous  new  forecasting  problems  ^Silence  dls- 

our  forecast  workers.  They  must  acquire  experience,  dls 

cover  basic  principles,  and  prepare  their 
relatively  efficient  techniques  have  gradua_ly  Ueen  in 
dustrilllLd."  Many  weather  forecasts  are  presented  num¬ 
erically  o 

In  the  future  the  relation  between  numerical  forecast- 

ing  Lrord?nary  dally  forecasting  Sro- 

relation  betweeti  mechanized  industry  and  handicraf  p 

ductions. 

At  present,  researches  in  numerical 
not  only  constitute  the  necessary  preparation  for  future 
numerical  weather  forecasts  but  also  considerably  aid  in 
anSyzlS  and  explaining  the  physical  causes  of  weather 

changes . 
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A 


with  this  In  mind,  we  initiated  dally  weather  forecasts 
during  the  development  of  the  people's  meteorological  en¬ 
terprises,  and,  in  1954,  we  started  researches  in  numerical 
weather  forecasting. 


II 


As  in  other  types  of  neteoi^ologlcal  work,  the  most 
basic  thing  Ih  numerical  itfsather  forecasting  is  to  under¬ 
stand  the  elements  of  the  atmospheric  process.  Such  an 
understanding  essentially  depends  on  the  analysis  and 
generalization  of  actual  atmospheric  changes.  The  neglect 
of  ,  these  elements  would  transform  numerical  weathei*  analysis 
Intb  metaphysics. 

Of  course,  the  area  of  our  work  is  not  confined  to  num¬ 
erical  weather  forecasting  alone.  We  have  completed  vital 
research  in  characteristics  of  the  atmospheric  process. 

Hsu  Erh-hao^  was  the  first  to  study  the  development  of 
the  East  Aslan  cold  front,  and  its  vorticlty,  extent  of 
divergence,  and  distribution.  In  addition,  he  also  cal¬ 
culated  the  value  of  each  term  in  the  humidity  equation. 

But  the  result  of  his  calculation  deviated  significantly 
from  calculations  made  by  the  usual  approaches. 

Researches  by  Peiping  University^ show  that,  under 
particular  conditions,  the  values  of  various  terms  of  the 
vorticlty  equation  are  similar  to  those  of  the  equilibrium 
equation.  Special  attention  must  be  directed  to  under¬ 
standing  these  deviations,  because  they  suggest  that  methods 
of  simplification  are  possible  for  certain  basic  equations, 
such  as  the  vorticlty  equation. 

As  far  as  the  existerice  of  the  eddy  transport  effect  in 
large-scale  general  circulation  is  concerned,  numerical 
analysis^  was  carried  out  on  the  basis  of  data  obtained  in 
February  1956.  The  res'.ilt  explained  various  terms  of  the 
humidity  equation,  in  regard  to  the  two-layer  geostrophlc 
model;  the  term  eddy  transport  is  very  broad  and  it  can¬ 
not  be  easily  Illustrated  by  the  mean  transport  field. 

Therefore,  it  is  Improper  for  the  capitalist  countries 
to  neglect  the  eddy  transport  term  or  to  replace  it  directly 
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by  the  mean  transport  field*  All  these  works  are  neces¬ 
sary  for  establishing  the  corresponding  models  of  numeri¬ 
cal  forecasting*  The  amount  of  this  kind  of  work  is  still 
limited,  and  further  development  along  this  line  is  neces¬ 
sary. 


III 


The  essential  prepaj'ation  for  numerical  weather  fore¬ 
casting  consists  of  the  followings  on  the  basis  of  past 
weather  experience,  we  apply  the  theory  of  dynamics  to 
study  some  simple  conditions  by  means  of  numerical  analysis, 
so  that  the  general  characteristics  of  the  effects  of  a 
given  physical  element  can  be  obtained.  For  example,  it 
is  commonly  known  that  an  eddy  produces  quantitatively 
various  transport  effects*  But  eddy  transport  can  reduce 
vortlcity,  a  phenomenon  similar  to  reduction  of  a  high 
field,  because  the  reduction  process  of  a  high  field 
naturally  causes  the  reduction  of  the  relative  vortlcity 
field 


It  proves  that  the  process' 


=  AV 


2yo 


(where  A  is  the  exchange  coefficient)  actually  causes  the 
reduction  process 

0*  :z  0  K 


and 


(1) 


where  At  and  As  are  the  time  and  apace  interval,  re¬ 
spectively.  Thus,  if  we  neglect  the  unnecessary  small  dis¬ 
turbances  and  the  accumulated  observation  errors,  cut-off 
errors,  and  other  probable  errors  which  add  to  the  smooth¬ 
ing  process  effect,  numerical  weather  forecasts  can  be 

*  Obykhov  did  similar  research* 
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achieved  vlth  the  Inherent  eddy  effect  of  the  model*  Tsao 
Chl-p*lns^  later  senerally  analyzed  these  properties  by 
means  of  the  stability  of  motion* 

Topography  may  greatly  affect  circulation*  It  plays 
an  especially  Important  role  in  China.  For  determining 
the  magnitude  of  its  effect  on  weather  development*  it  Is 
necessary  to  make  cert eih  calculations  on  the  basis  of 
assumed  simple  and  ideal  conditions •  Under  conditions  of 
general  west  wind  ahd  north  wind,  graphical  methods®  were 
used  to  obtain  the  result  of  topographical  effects  on 
general  baro tropic  oiz*culatl6n  disturbances* 

The  numerical  tendency  results  obtained  ind lent e  that 
the  topographical  disturband^  is  mubh  greater  then  that  of 
liionth  America  under  similar  conditions*  Under  baroclinlc 
pressure*  because  of  the  effect  of  the  stable  layer*  the 
topographical  disturbance  of  the  lower  layer  is  greater* 

In  the  mountain  region  of  A-erh-t'al  Shan  and  Hsiang- 
t*len  Shan'*  if  the  mean  west  wind  velocity  is  five  meters 
per  second*  the  high  t|ndenoy  at  700mb  on  the  anterior 
mountain  slope  may  exceed  40  meters  in  12  hours*  Besides* 
with  barotropio  models,  calculations  were  made  on  the  ten¬ 
dency  of  topographical  disturbance  in  north  China  under 
general  north-west  circulation.®  It  was  discovered  that 
the  position  and  intensity  of  Y1  Chao  in  north  China  cor¬ 
responded  to  the  initial  stage  of  north-west  circulation* 
All  these  activities  provide  us  with  a  conception  of  topo¬ 
graphical  disturbances  under  bauratroplc  and  baroclinlc 
pressures  in  China* 

Under  actual  baroclinlc  pressure*  the  basic  force  is 
sustained  and  supplied  by'outslde  heat  sources*  Therefore* 
for  medium  and  long-range  weather  process  forecasting* 
the  process  of  how  the  heat  sources  affect  the  atmosphere 
is  an  important  problem. 

In  1957,  Tsao  Chi-p*lng^  and  others  used  the  linear 
two-layer  model  to  calculate  the  effect  of  heat  sources 
under  ideal  conditions*  They  discovered  that*  no  matter 
whether  the  heat  source  was  considered  as  constant  or  In¬ 
fluenced  by  the  temperature  and  pressure  field  (reverse 
process),  the  effect  of  heat  sources  could  be  adjusted 
after  five  days  or  two  weeks*  implying  that  the  temperature 
and  pressure  field  lagged  behind  by  about  half  a  month* 
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But  because  of  difficulties  ia  designins  the  models,  the 
effect  of  calculation  with  models  is  Insufficient. 

Chu  Pao-chen  and  others^®  tried  the  Kibel  radiation 
method.  They  introduced  thh  tiictCr  6f  had!  at  ion  into  the 
linear  two-layer  models  and  fihhlly  obtained  the  equation 
of  the  heat  source  eff^cst  on  the  high  field.  Actual  cal¬ 
culations  were  carried  out.  Of  great  significance  la  the 
fact  that,  when  considering  the  relation  (exchange  coeffi¬ 
cient)  between  Input  (heat  source)  and  output  (high  field), 
the  state  of  the  heat  source  when  affected  by  the  tempera¬ 
ture  and  pressure  field  Is  the  same  as  the  state  when  the 
heat  source  is  constant. 

The  final  basic  forms  of  the  different  influences  are 
the  same.  The  calculation  results  of  both  situations  are 
alike.  In  addition,  Chu  Pao-chen  and  others  studied  the 
instability  of  the  sine  curve  and  discovered  that  in  such 
a  situation  the  instability  should  be  determined  by  more 
abundant  data.  All  these  researches  provide  us  with  a  con¬ 
ception  of  the  effect  and  the  actual  process  when  a  heat 
source  is  applied  to  a  temperature  and  pressure  field.  At 
the  same  time,  it  creates  problems  of  model  design. 


IV 


The  realization  of  numerical  forecasting  depends  upon 
actual  calculation  with  models.  Evidently  a  great  amount 
of  work  must  be  done  in  establishing  and  choosing  models, 
and  in  developing  methods  to  obtain  solutions.  Only  by 
means  of  actual  calculations— the  process  of  verifying 
theory  by  fact,  may  we  reach  a  better  understanding  of 
forecasting  by  models. 

As  far  as  the  one-layer  quasi-geostrophlc  model  of 
large-scale  motion  is  concerned,  graphical  methods  can 
be  used  in  actual  forecasting.  For  this  reason,  at  the 
initial  stage  (1954-1956),  we  carried  out  some  researches 
in  graphical  numerical  forecasting. 


Liao  Tung-hsien^^  examined  the  two-layer  forecast 
equation 

|S  *  _  (2) 

where  ^  denotes  absdlute  yortiQlty»  'Vk  denotes  thermal 
wind,  and>^^r  denoted  th#  mSan  flow  v’drticlty  of  the  thermal 
wind*  Oalculatlon  proves  that  tihe  second  tehm  on  the  right- 
hand  side  of  tlie  equation  IS  k  i‘eldtiVeiy  smkll  bei^rectlon 
factor*; 

In  primary  approximation*  it  is  reasonable  to  hSSUme 
that  the  correction  factor  varies  insignificantly  duf'itig 
the  forecast  period*  On  this  basis  the  approximate  solu** 
tion  can  be  obtained  graphically  according  to  'V'  and  Vr 
by  manipulating  the.^  and /77-  field  separately*  Actual 
forecasts  show  that  better  results  are  attained  when  we 
apply  this  method  to  forecasts  of  cold  front  development 
than  [when  we  apply  it]  to  forecasts  of  the  mean  flow  of 
nonthermal  driven  wind  (barotroplc  model)*  This  method 
is  basically  similar  to  that  of  Pjortofts. 

1  p 

Subsequently*  a  researcher,  by  adapting  a  similar 
method, .tested  a  48-hour  500mb  Isobaric  surface  forecast. 

It  shows  that  some  cold  fronts  can  be  forecast  48  hours  in 
advance.  However,  the  blocking  high  pressure  cold  front 
cannot  be  forecasted  with  certainty.  Besides,  since  the 

y§§. 

term  is  relatively  small,  It  can  be  neglected  if  graphical 
solutions  are  applied  in  numerical  forecasting* 

Bolin  used  the  two-layer  model  graphical  method  to 
obtain  the  solutions  of  the  situations  of  both  layers* 

Hsu  Erh-hao^5  applied  this  method  to  the  Sawyer  and  Bushby 
model  and  obtained  the  corresponding  forecast  equation, 
which  is  being  verified  experimentally**  Several  experi¬ 
ments  show  that  the  results  are  superior  to  those  obtained 

♦Recently,  Hsu  Erh-hao  established  a  l,OOOmb  direct 
graphical  forecast  method  to  be  published  in  the  near 
future. 
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by  using  the  Ellasaen  and  Reed  Oodel.  But  the  work  in¬ 
volved  Is  more  complloated* 

In  two-layer  models »  a  propoi^tibhni  deenease  dt  temper¬ 
ature  generally  Is  oonst^inti  iSln^ij  thd  upt)er  and  lower 
layers  of  the  temperaturf  iTlhid  are  similar.  It  Is  Impos¬ 
sible  to  consider  the  existence  of  the  front  poslton.  Be¬ 
cause  of  this  fact,  application  of  the  general  atmospheric 
pressure  coordinate  does  not  yield  good  results  In  fore¬ 
casting  weather. 

Gh*en  Hslung-shan^^  and  others  pointed  out  that  the 
existence  of  the  front  position  can  be  considered.  If  the 
two-layer  model  of  the  position-temperature  •  coordinate 
Is  used.  Because  the  distance  of  the  front  position  per 
equivalent  surface  is  particularly  small,  wo  can  consider 
the  equal  d  areas  above  and  below  the  front  position  as  the 
required  layers  and  apply  the  upper  and  lower  limits  to 
it;  In  this  way  we  are  able  to  Indicate  the  difference  of 
stability  both  within  and  without  the  front  poslton.  The 
flow  field  of  the  9  surface  can  thus  be  forecasted.  They 
considered  the  upper  limit  ©3  s  constant  and  the  lower 
limit  ©2  =  constant.  The  positional  vortlclty  s 

""  1,566'  ' 


where  K  =  (Op  -  C9)/Cp,  denotes  the  flow  function  of  9 
surface,  and  others  are  conventional  symbols. 

Based  upon  the  above  method,  a  graphical  forecast  was 
made  (12  hours).  The  result  was  found  to  be  relatively 
better.  Xet  the  difflciilty  of  this  method  lies  in  the 
problem  of  assuming  the  boundary  conditions,  especially  the 
boundary  condition  of  the  upper  limit.  Proper  solutions 
of  this  problem  remain  to  be  found.  Besides,  the  existence 
of  these  front  positions  has  not  been  considered  in  con¬ 
structing  models* 

Besides  the  experiments  In  the  graphical  method,  the 
accuracy  of  Piortoft*s  graphical  method  Is  being  studied. 
Pjortoft  considered  the  difference 

(3) 
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In  Poisson’s  equation 

^  _ 

where  ^  ,  Cf^  is  the  mean 

value  of  (X  Inside  the  reotanglis  and  d  denotes  the  inter¬ 
val  of  the  rectangle*  Then  thd  solution  ia 

C?<  =  (1  +  M  +  .•.«...*)h.  (4) 


It  Is  aOcurate  only  if  an  infinitive  region  la  assumed . 
Fjortofi  proved  that 

CX  ^  =  h  +  2fi* 


can  be  considered  as  the  first  approximate  solution  of 

(3). 

Ch’ou  Ghi-fan^^  and  others,  in  thirty  experiments, 
proved  that  error  in  using 

CtC  =  d  +  2h 

is  within  the  permissible  degree  of  accuracy  of  the  gra¬ 
phical  method.  But  h  consists  of  the  vortlclty  mean  cir¬ 
culation.  It  is  obtained  by  graphically  extrapolating  the 
mean  circulation  in  space.  But  based  upon  the  same  ex¬ 
periment,  with  the  distance  of  the  mean  circulation  rec¬ 
tangle  d  =  600  millimeters,  the  change  of  value  is  as 
large  as  two  thirds  the  original  circulation  field.  The 
result  is  therefore  unsatisfactory. 

The  graphical  method  for  numerical  forecasts  of  rain 
is  also  necessary.  Chang  iCen-*-'  suggested  a  synthetic 
process  for  large-scale  short-range  numerical  forecasts 
of  rain.  Particular  examples  were  chosen  for  purposes  of 
calculation,  and  the  results  wore  good.  The  method  la 
worthy  of  further  investigation  and  Improvement, 

For  accurate  numerical  forecasts,  the  use  of  eleoti*onlc 
computers  is  necessary?  hence,  the  finite  differential 
method  is  used.  In  this  connection,  one  should  transform 


8 


the  differential  equations  of  numerical  weathet*  forecasts 
Into  differential  equations  for  further  calculation*  A 
lars®  number  of  techniques  are  used  in  all  these  numerical 
Calculations,  whiph  InvolV'e  the  width  6f  the  rectangle, 
the  interval  of  tim^,  and  the  brbblem  bf  SXtrabolatlon  of 
time.  Since  no  actual  caipuiatlbil  has  been  performed  by 
electronic  computers,  no  final  conclusion  has  been  drawn, 
as  yet,  in  this  respect. 

As  far  as^golsson*s  approximation  method  is  concerned i 
Ku  Chen-tsao^®  used  the  least  square  method  and  a  finite 
number  of  the  mean  value  terms  M  to  obtain  the  approximate 
solution  of  equation  (3).  He  proved  that,  when  choosing 
h  and  h,  the  best  approximation  form  is 


CA*  -  0.4h  +  4,4  h,  4  ^  2/d  ^  12;  (6) 


This  is  more  accurate  than  Pjortoft*3  formc?^*  s  h  +  2h  or 

=  h  +  3  h.  It  is  also  simpler  and  more  accurate 
than  Belousov *s 

=  (3  h  +  2h  +  E),  (7) 

It  is  convenient  to  use  approximate  formula  (6)  in  elec¬ 
tronic  calculations.  For  a  higher  degree  of  precision# 
similar  methods  can  be  used  to  obtain  that  mean  value  of 
h  which  may  yield  the  best  approximation  formula. 

The  degree  of  precision  of  the  graphical  method  is 
limited.  In  general,  it  is  much  more  accurate  to  use  the 
numerical  approximation  solution  of  the  vorticity  equation. 
But  the  application  of  different  methods  of  numerical  ap¬ 
proximation  involves  problems  of  the  amount  of  work,  the 
degree  of  accuracy,  and  the  stability. 

Therefore,  it  is  necessary  to  study  further  the  pro¬ 
blems  of  obtaining  solutions^ to  the  vorticity  equation. 

A  Mathematics  Institute  team^^  found  the  circulation  co-^ 
efficient  of  Bllnova*8  barotroplo  linear  forecast  equa¬ 
tion  as  a  logarithmic  solution  In  terms  of  time  t  by  trans¬ 
forming  the  original  differential  equation  into  an  integral 
equation  under  conditions  of  a  longitudinal  function. 
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Another  method  of  solutlon^^  has  also  been  considered. 

It  Is  worthy  of  sehlous  consideration  and  experimental 
verification* 

Like  the  graphical  method*  when  manipulating  extra¬ 
polation*  the  ohe^lkyer  model  is  "the  most  simple  to  apply. 
Based  on  two  twelve-hcur  experimental  forecasts,  Wang 
Yao-sheng  and  othCrs^l  fouhd  that  during  the  summer  In 
China  It  Is  better  to  use  ihe  bne-layer  model  for  TOOmb 
than  for  500®b, 

However i  the  one*i'iayer  model  is  always  insufficient. 

For  more  genehal  application,  Liao  Tung-hsiett  suggests 
a  two-layer  model  to  be  applied  to  situations  of  parti¬ 
cular  distribution  of  vertical  motion.  He  divided  the 
atmosphere  Into  two  layers,  which  were  described  separate¬ 
ly  by  a  two-layer  model.  Prom  this  two  -  layer  model, 
the  three-layer  model  forecast  equations  were  derived. 

It  la  convenient  to  obtain  the  solutions  graphically 
from  these  equations,  A  weather  bureau  team  found  an 
Infinite  aeries  solution  for  the  two-layer  model  fore¬ 
cast  equation  expressed  In  terms  of  spherical  coordinates. 

In  certain  sltimtlons,  the  three-layer  models  are 
essential.  In  various  regions,  If  the  temperature 
changes  of  the  upper  and  lower  layers  do  not  assume  a 
linearly  proportional  development,  the  use  of  two-layer 
models  for  weather  forecasting  generally  Is  Impossible, 

In  fact,  the  blocking  high  pressure  developed  In  the 
Ural  region  produces  the  more  complicated  barocllnlc 
effect^?.  Thus,  three-layer  models  are  required  for 
forecasting. 

In  addition,  some  characteristic  rules  have  been 
suggested  in  aplylng  three-parameter  modela^^.  For 
excuzipl^f  th6  sp0ol3.1  factors  of  the  dovsloptaont 
three-layer  models  can  be  observed  from  the  proportionate 
decrease  rate  of  the  mean  circulation  of  thermal  wind. 

If  high  altitude  multiple-layer  weather  charts  are  pro¬ 
vided,  the  following  rule  can  be  applied  to  forecast 
qualitatively  the  barocllnlc  development,  which  Is  a 
particular  property  of  the  three-layer  models^.?:  In  an 
extensive  region  the  cold  (warm)  mean  circulation  In¬ 
creases  (decreases)  In  strength  In  an  upward  direction; 
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hence,  the  altitude  of  the  Intermediate  iisobarlc  layer  Is 
elevated,  or  otherwise  Ibwered. 

Simplified  eb^lmation  pi^buiatibhs  show  that.  In  China, 
topographic  effiots  iiuBt  ob  considered  when  using  numeri¬ 
cal  forecast  modbis.  Nevertheless,  how  to  describe  exactly, 
and  illustrate  by  models,  the  topographic  effects  remains 
a  highly  important  problem* 

Theoretically,  GreenS  function  can  be  used  to  solve 
the  problems  of  barocllnic  numerical  weather  forecasting 
under  topographic  Influence*  Unlike  Sh*  A#  Musayelyan, 
who  considered  sea  level  as  the  boundary  condition,  Wang 
Chung-hao^^  considered  the  local  region  as  the  boundary 
condition  In  solving  the  forecast  equation.  The  results 
show  that  Green's  functions  (Including  the  Integrals  of 
Green's  functions)  thus  obtained  are  regionally  deter¬ 
mined. 

Hence,  the  value  of  Green's  function  obtained  at  one 
point  Is  different  from  that  obtained  at  another  point. 
Though  geostrophic  wind  Is  assumed,  the  solutions  of  the 
equations  are  highly  complicated.  It  Is  obvious  more 
difficult  to  obtain  a  complete  solution. 

Lin  Hslng-yuan^^  used  similar  boundary  conditions  and 
suggested  the  Fourier  Method  for  obtaining  solutions  of 
the  forecast  equations.  He  assumed  the  topographic  boun¬ 
dary  condition 

=  1  -Af, 

where  Is  proportional  to  the  geometric  topographic 
altitude. 


{^=  P/Po,p,  , 


p  denotes  surface  atmospheric  pressure,  Po  denotes  stan¬ 
dard  sea  level  atmospheric  pressure }•  Instead  of  using 
the  parameter,  he  expanded  the  solution  Into  logarithmic 
series  In  X  »  found  the  J-th  order  logarithmic  forecast 
equation  In  A  »  the  j-th  order  approximation,  and  thus 
obtained  the  method  of  Its  solution.  He  denoted  the  ver¬ 
tical  coordinate  of  the  equation  by 

^  »  p/p 
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(where  P  Is  surface  atmospheric  pressure,  p  is  high  alti¬ 
tude  atmospheric  pre^aure),  so  that  the  lower  limit 

^  =  1» 

He  applied  hexagbha  io  otitaln  the  approximate  Colutlon  of 
the  fCil^edaBt  equatlCh*  l?he  daloulation  also  considers 
topographic  effects.  According  to  ijiore  than  40  examples 
applied  lii  hdhth  China,  a  consideration  of  topographic 
effects  has  raised  the  degree  of  accuracy  by  ten  percent# 

Under  barocllnlc  conditions,  the  problems  of  calcula¬ 
ting  topographic  tendencies  have  increased  in  number • 
actual  tendency  calculations  show  that  the  accurate  deter** 
mlnation  of  barocllnlc  topographic  disturbance  tendency 
is  much  more  complicated  than  previously  imagined.  The 
first  problem  is  how  to  determine  the  surface  slope  wind. 

If  we  simply  consider  850mb  or  the  upper  limit  bound¬ 
ary  layer  wind  as  the  basis  for  calculating  the  vertical 
motion  of  the  elope  wind,  we  may  not  be  able  to  determine 
the  expected  topographical  tendency.  Actually  the  rela¬ 
tionship  between  surface  wind  and  lower  layer  wind  of  the 
free  atmosphere  is  not  so  simple*  If  we  do  not  depend  on 
actual  observations,  we  find  it  difficult  to  determine  the 
surface  wind  and  its  time  rate  change  (wind  ’’tendency  ) 
by  using  the  equation  of  motion.  If  the  degree  of  accur¬ 
acy  cannot  be  determined  by  actual  observations,  the  pro¬ 
per  choice  of  boundary  conditions  for  forecast  calcula¬ 
tions  may  become  a  problem.  Thus  it  la  necessary  to  per¬ 
form  the  calculations  more  rigorously,  to  gain  more  ex¬ 
perience,  and  to  delve  more  deeply  into  theoretical  an¬ 
alysis. 


V 


In  early  1958,  with  the  earnest  and  unselfish  assis¬ 
tance  of  Soviet  specialist  Dobryshman,  who  was  motivated 
by  internationalism,  the  Central  Weather  Bureau  initiated 
medium  and  long-range  numerical  forecasting.  They  did 
much  of  the  preparatory  work,*^°*  ^  and  attempted  long- 

range  forecasts.  Keeping  pace  with  the  leap  forward 
movement  of  1958,  numerical  forecasting  took  many  bold 
steps. 
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Various  work  groups  suggested  many  medium  atxd  long- 
range  linear  for^eoadt  mddels  and  tasted  certain  oaloula- 
tlons.  Thld  Is  ajddrrect  l&eglnnlugi  because  for  medium 
and  long-range  fOfeekstiJ  assisiatice  of  new  methods  Is 
urgently  needed.  It  Is  imjpdfsl^le  to  determine  the  merits 
and  defects  of  all  models  wishoUt  bold  experiments  and 
practical  verifications. 

Chu  Pao-ohen  suggested  a  two-layer  geostrophlo  model: 
constant  heat  source,  as  well  as  a  geostrophlc  model  which 
takes  into  cohsideratlon  radlaktion  and  vortex  heat  (deal^ 
Ing  with  Inversion  caused  by  the  temperature  pressure 
field }•  The  Klbel*  method  is  followed  In  considering  the 
radiation  aspect  (referring  to  the  corresponding  non¬ 
linear  three-dimensional  geostrophlc  model Certain 
properties  of  these  models  were  understood  after  some 
linear  calculations  were  carried  out. 

Besides,  five-day  forecasts  and  one-month  forecasts 
were  tested  by  using  the  barotroplc  topographic  mode. 
Useful  experiences  have  been  obtained.  All  these  works 
have  helped  solve  the  mystery  of  long-range  numerical 
weather  forecasting,  and,  by  their  means,  experiences  In 
calculation  were  obtained. 

Our  large-scale  calculations  of  linear  models,  the  re¬ 
sults  of  experiments  of  calculations  of  linear  models 
carried  out  by  foreign  countries  (Blinova,  Mashkovlch), 
and  the  problems  educed  during  these  calculations  natur¬ 
ally  lead  us  to  consider  under  what  conditions  the  use  of 
linear  models  are  permissible  and  to  consider  the  methods 
of  designing  the  models  and  determining  the  parameters 
under  these  permissible  conditions,  so  that  the  stability 
and  superiority  of  the  models  can  be  guaranteed. 

Researches  have  been  conducted  for  this  purpose.  Re¬ 
cently  new  appiToaches  have  been  adopted  to  analyze  further 
the  problem  of  stability  of  models,.  Relatively  more 
general  results  have  been  obtained,^  This  greatly  aids 
In  designing  linear  models. 

In  short,  this  kind  of  work  is  at  Its  Initial  stage. 

Kie  possibilities  of  exploration  are  unlimited.  Rese¬ 
arches  of  similar  nature  can  be  conducted  In  dealing 
with  non  linear  problems. 
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VI 


Theoretically  speaking*  the  more  basic  problem  is  not 
the  determination  of  models i  but  the  method  of  approach 
used  In  nunlerloal  forecastsi  Since  V,  Bjerknes,  numberl- 
cal  weather  forecasting  has  been  considered  to  be  ® 
blem  of  mathematic problem  of  obtaining  the  initial 
values  for  solving  the  system  of  partial  differential 
equations*  ‘ 

That  is,  one  should  use  the  initial  value 
partial  derivative  of  the  meteorolbglcal  quantity  to  find 
the  solution  Of  the  system  of  differential  equation*  But 
dally  Weather  forecasting  is  based  on  the  historical 
development"  of  the  change  of  weather  (namely  the  recent 
period,  within  several  hours  or  several  days). 


the  long-range  weather  forecasting 
these  two  approaches  are  completely 
have  their  own  particular  merits, 
forecast  results.  The  problem  now 
any  correlation  between  these  two 


It  is  similar  to 
approach.  Although 
different,  but  they 
Judging  from  actual 
is  whether  there  is 
approaches. 

It  is  obviously  improper  to  think  of  dally  weather 
forecasting  as  simply  a  Meduotive  process,  and  not  to 
ftlve  it  serious  consideration.  As  a  matter  of 
under  certain  conditions,  it  is  possible  to  prove  that 
these  two  appraches  are  identical.  For  large-scale  cir¬ 
culation  under  "approximation  of 

conditions,  Ku  Ohen-tsao  has  proved  that:  in  considering 
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(II) 


(where  Is  the  altitude  of-  equal  pressure  surface,  J 
la  the  Jacobian, 

r  at 


Is  the  altitude  tendency,  cK  Is  a  known  parameter,  V 
denotes  pressure,  and  f  Is  geostraphlc  parameter); 
.the  condition 


♦(»,  y,  y,  *)  r»  #*) 

for  all  t  -'"® 


(12) 

(IJ) 


(where  t^  denotes  Initial  time,  P  denotes  a  known  func¬ 
tion)  and  condition  when  P  = 


when 


♦(*,  y,  0  “  ^i(*i  y»  <)»l  I 

2*(*,  y»  M  0  y»  <)** 


(14) 


! 

I  =  known  distribution,  (15) 


and  when  <•>/  —  «  ,  p»0,  •■'0  (i«^ 

a 

(where  t  -Cf-<  t^)  are  equivalent  In  value.  In  (14)  P^ 
can  be  replaced  ny  other  values  of  pressure  P  ss  P^. 
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For  weather  forecasts  within  limited  areas »  because  of 
the  effect  of  sea  level  boundary  conditions,  it  Is  neces¬ 
sary  to  revise  slightly  the  equivalent  donditions  al¬ 
though  they  are  basically  similar* 

Tide  above-mantiohed  equivalent  value  property  shows 
that  the  change  of  the  suiffaoe  temperature  pressure  field 
suggests  the  three-dimensional  temperature  pressure  field 
structure  of  the  bat^oclihic  atmbsphere  and  determines  the 
development  of  the  three-dimensional  temperature  pressure 
field.  Because  of  the  baroolinlc  property,  we  have  to 
ootisider  the  “historical  development"  of  weather  when  us¬ 
ing  the  second  approach* 

But  the  equivalent  value  problem  Is  only  a  particular 
problem.  Actually  these  two  approaches  do  nqt  yield  equi¬ 
valent  values*  To  state  this  explicitly,  historical  data 
are  Indispensable  for  any  numerical  weather  forecast*  so 
that  errors  may  be  reduced* 

By  merely  observing  a  few  differential  equations  we 
can  not  fully  understand  this  problem.  Numerical  weather 
forecasting  Is  not  purely  a  mathemetlcal  problem.  Prom 
our  experiences  we  understand  the  necessity  of  "historical 
data,"  Like  the  principle  of  rhythm  we  have  experienced, 
the  conditions  of  various  periods  In  the  past  play  a  con¬ 
siderably  Important  role  In  long-range  forecasts. 

Of  course,  as  far  as  the  general  forecast  problems 
are  concerned,  the  most  Important  problem  Is  the  Inade¬ 
quacy  of  observational  procedure;  hence,  we  are  unable 
to  avoid  errors.  Past  data  are  abundant* 

Therefore,  the  problem  of  choosing  the  correct  approach 
should  be  based  upon  the  availability  of  observation 
data;  we  should  choose  an  approach  so  that  available  ob¬ 
servational  data  can  be  fully  utilized,  and  we  should 
avoid  approaches  which  require  unavailable  observational 
data.  For  this  reason  It  Is  advisable  to  utilize  data  In 
the  past,  and  as  a  matter  of  fact  we  cannot  avoid  using 
them. 

However,  further  research  Is  necessary  on  how  to 
utilize  fully  the  most  reliable  and  complete  distribution 
of  meteorological  data  in  past  and  recent  periods,  so 
that  the  most  proper  form  of  the  system  of  weather  equa- 
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t ions,  when  Integratihg  numerically  with  respect  to  time, 
can  be  chosen  and  the  proper  oonditiohs  foi*  the  solution 
can  be  diitermin^*  Researches  in  thic  t'ecpect  are  at 
their  initial  sthgei  The  fiblutibn  of  this  problem  may 
provide  a  new  direction  fot^  numerical  weather  forecasting* 


# 


* 


* 


Certainly,  the  history  of  the  development  of  numerical 
weather  forecasting  in  China  is  short*  However,  under 
the  correct  leadership  of  the  Party  and  the  State,  the 
development  of  numerical  weather  forecasting  in  such  a 
brief  period  has  been  comparatively  rapid.  Actually, 
only  after  the  establishment  of  the  people  s  meteorological 
enterprises , and  of  the  complete  observation  and  forecast 
system,  were  the  foundation  and  necessary  conditions  for 
all  the  researchers  provided. 

During  the  period  of  reactionary  rule,  all  these  c^- 
dltlons  were  incomprehensible.  Since  the  leap  forward 
movement,”  numerical  weather  forecasting  in  China  has 
assumed  a  new  appearance.  It  is  expected  that  in  the 
near  future,  by  combining  practical  requirements  with 
new  technology  and  new  theory,  numerical  weather  fore¬ 
casting  will  reach  a  higher  stage. 
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